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Abstract 

The electric field increases toward infinity in the narrow region between closely 
adjacent perfect conductors as they approach each other. Much attention has been de- 
voted to the blow-up estimate, especially in two dimensions, for the practical relevance 
to high stress concentration in fiber-reinforced elastic composites. In this paper, we 
establish optimal estimates for the electric field associated with the distance between 
two spherical conductors in n— dimensional spaces for n > 2. The novelty of these 
estimates is that they explicitly describe the dependency of the blow-up rate on the 
geometric parameters: the radii of the conductors. 

MSC-class: 15A15, 15A09, 15A23 



1 Introduction 

We consider the blow-up of the electric fields in the narrow region between a pair of per- 
fect conductors which is closely adjacent in n dimensions (n > 2). Conductors provide 
higher intensity of electric flux around them. The intensity increases as a pair of con- 
ductors approaches each other, and the electric field even reaches toward infinity (refer to 

In this paper, we present the optimal blow-up estimate for the electric field with respect 
to the distance between a pair of conductors under the assumption that the conductors are 
of spherical shape in n dimensions (n > 2). The novelty of these estimate is to describe 
explicitly the dependency of the blow-up rate on the radii of the conductors: this paper is 
the first result to establish the role of the geometrical factor of conductors in the blow-up 
of the electric field in three or higher dimensions. 

Besides the consideration of the gradient estimates in the frame of the electrostatic the- 
ory, much attention has been drawn to it of the relevance to the stress-strain behavior of 
composite materials, especially in two dimensions. According to Budiansky and Carrier 
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[5], unexpectedly low strengths in stiff fiber-reinforced composites have been reported, due 
to the high stress concentration occurring in the narrow region between fibers (also refer 
to [5]). In the anti-plane shear model, the stress tensor represents the electric field in the 
two dimensional conductivity model, where the out-of-plane elastic displacement satisfies 
a conductivity equation [5J. Thus, the gradient estimates for electric field have a valuable 
meaning in relation to in the failure analysis of composite material. To give a brief descrip- 
tion of related works, for the case that the inclusions and the outside of inclusions have 
the comparable conductivities (or shear moduli), it was verified that the electric field re- 
mains bounded independently of the distance between the inclusions. Li and Vogelius [TU] 
have shown that the electric field does not blow up even when the inclusions approaches 
each other. Moreover, Li and Nirenberg [5] have extended this result to elliptic systems. 
These results point out that the extremely high conductivity (or the stiffness of fibers) is 
indispensable to the blow-up phenomena. 

In this respect, much attention has been focused on the model of a pair of perfect 
conductors which are e apart. Ammari, Kang, H. Lee, J. Lee and Lim [21 [3] have established 
the optimal blow-up rate e -1 / 2 as the distance e goes to zero, when conductors are of 
circular shape in two dimensions. Yun [TTJ has extended the above mentioned result 
to a sufficiently general class of the conductors' shapes in two dimensions. In three or 
higher dimensional case, Bao, Li and Yin [1] recently obtained the optimal blow-up rate 
for perfect conductors of general shape: the optimal blow-up rate is (e|loge|) _1 for three 
dimensions, and is e _1 for higher n dimensions (n > 4). However, their estimates are only 
given by the distance between two conductors and geometric information of conductors are 
not incorporated into the estimates. 

What is new in this paper is that for the case of spherical perfect conductors in three 
and higher dimensions, the gradient estimates are established in terms of the radii as well 
as the distance between inclusions. What is more is that the approach introduced in this 
paper to derive the estimates is distinct from the methods of Bao et al. [I] and Ammari 
et al. In the two dimensional case, our approach provides the same estimates as of 

Ammari et el., Proposition 13. 21 in a much simpler way for the case of perfect inclusions. 



2 Mathematical formalism and main results 

From now on, R™ denotes n dimensions, and B r (x\, X2, ■ • • , x n ) is the sphere with radius r 
and center (xi, X2, • • • , x n ) in R™. Given any entire harmonic function H in R™ (n > 2), we 
define the electric potential u as the unique solution to the following conductivity problem: 

' Au = 0, in R"\(Di UD 2 ), 

< U (x)-i/(x)=0(|x| 1 -") as|x|^oo, 
u\dDi — Ci (constant), 

IdD ® vU = 0, f° r * = L 2, 

where x = (x\,X2, • ■ ■ ,x n ). This solution u can be interpreted physically as the electric 
potential outside conductors D\ and D2 under the action of applied electric field VH . 

In this paper, we start by considering the case that V H is a uniform field, i.e, H = a • x 
for some constant a in R™, in Theorem 12.11 and 12.31 Based on these, the optimal upper 
bound of the gradient for any entire harmonic function H is established in Theorem l2.4l 
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Theorem 2.1 (Three dimensions) We assume that Di and D 2 are the pair of spheres 
with radii T\ and r 2 that are 2e apart in R 3 . Thus, we set 

D 1 = B ri (n + e, 0, 0) and D 2 = B r , 2 (-(r 2 + e), 0, 0). 

Let u be the solution to for H(x\,x 2 ,x^) = ^2 i=1 CLiXi- Then, there exists a positive 
constant C* independent of e, n, 7*2 and (01,02,03) such that 



u\dD! ~ u\ dD2 



<akx|(^p-]|loge| 
n + r 2 



for sufficiently small e > 



(a) In the case that a\ is nonzero, for any sufficiently small e, there is a point xo between 
D\ and D 2 such that 

^r\ax\ (J^L-) -J— < |V«(xo)|. 
2C* \ r i + r 2/ |e log e| 

The lower bound above is optimal in the sense that there is a positive constant C* 
independent of e, n, r 2 and (01,02,03), satisfying that 

I f\ 7*2 \ 1 

\\Vu\\ L ~ {R3 \ {DlUD2)) < C*\ai\ [—^ J po^' ( 2 ) 

for sufficiently small e > 0. 

(b) In the case that a\ is zero, the gradient of u does not blow up even when the distance 
e goes to zero, i.e., there is a positive constant Cq independent of e and (01,02,03), 
satisfying that 

\\Vu\\L°°(tt3\(D 1 UD 2 )) < Cq (|o 2 | + |o 3 |) , (3) 
for sufficiently small e > . 

Remark 2.2 The constant Cq at ([3]) depends on r\ and r 2 : in details, there is a constant 
C so that 

{r\ 7*2 
— , — 

when ai = 0. The derivation of the inequality above is included in the proof of Theorem \2.1\ 
The term 'Cmax|^-, ' becomes arbitrarily large for small r\ or 7*2, and then it is not 

guaranteed that the bound ^ above is optimal. However, our attention is focused on the 
contribution of r\ and 7*2 to the blow-up rate associated to the distance It. As mentioned 
in Theorem \2.1l the gradient of u blows up as e — > if and only if a\ is nonzero. On this 
occasion, the l/|e log e| term at ([2]) dominates the upper bound of ||V'tt||£«'(R3\(z>iU£>a))j an d 
the bound ((21) describes the contribution of r\ and 7*2 well. Therefore, Theorem \2.1\ would 
suffices for our purpose. 

Theorem 2.3 (Higher dimensions) We assume that D\ and D 2 are the pair of spheres 
with radii r% and r 2 that are 2e apart in IR ra (n > 4) . Thus, we set 

Di =B ri (n + e,0,--- ,0) andD 2 =B r2 (-(r 2 + e),0,--- ,0). 



(M + M), (4) 
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Let u be the solution to (QJ) for H(x\, X2, • • • , x n ) = Xa=i a i x %- Then, there exists a constant 
C** independent of e, n, r 2 and (ai,o 2 , • • • ,a„) smc/i that 



C** V r l+ r 2 



<C,Joi| 



7*1 ?" 2 



for sufficiently small e > 



(a) In the case that a± is nonzero, for any sufficiently small e, there is a point Xo between 
D\ and D2 such that 

Kl (^r) ~ < |V«(xo)|. 



2C„ 1 \ r i + 7 *2 / e 

TTie lower bound above is optimal in the sense that there is a positive constant C* 
independent of e, n, r-i and (01,02, ■ ■ ■ , a n ), satisfying that 



I|Vm|| L oc( R 3\( DiU£ , 2)) < C**|oi| ^ 



rir 2 \ 1 



T\ + T2 J e 



for sufficiently small e > 0. 



(b) In the case that a\ is zero, the gradient of u does not blow up even when the distance 
e goes to zero, i.e., there is a positive constant Cq* independent of e, r\, ri and 
(ai,a 2 , ■ ■ • ,a n ), satisfying 

n 

I|Vm|| L =c (R 3\ (£ , iU£)2)) < Cq* ^ |Oj|, 

i=2 

for sufficiently small e > . 

Similarly to Remark 1 2. 2 1 the constant Cq* above depends on n and r2,i,e, there is a constant 
C so that 



{r r 1 
r2 Tl ' i=2 



when ai = 0. The derivation of the inequality above is also presented in the proof of 
Theorem O 

In this work, the blow-up estimates in terms of radii is presented only for three or higher 
dimensional case. Speaking of two dimensions, Ammari et al [5] already provided the optimal 
bound ([5]) and (JTUJ) in terms of radii of circular inclusions as 



maxMu-H)\<C*l-^ — \= (5) 
w n + r 2 Je 



where r\ and r 2 are the radii of circular inclusions. This is also derived in Proposition 13.21 
of this paper. As has been mentioned before, the method in this paper is much simpler 
method. 
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Theorem 2.4 (General entire harmonic function H) Let Di and D2 be a pair of balls 
as assumed in the previous theorems in n dimensions (n > 3). We choose a large bounded 
domain Q containing D± and D 2 for any small e > 0. For the sake of convenience, we select 
the ball B^ ri+r2 ^(0, ■ • • , 0) as Q. For any given entire harmonic function H, let u be the 
solution to (QP for H . 

In three dimensions (n — 3), there is a constant C* , independent of e, r\, r2 and 
(ai,a 2 ,a 3 ), satisfying 



\Mu-H)\\ L ^ M3x{DlUD2)) < C*||V17|| £ ~ (n) (^^) pj; 
and 

( T\T2 \ 1 

Vu| \L<*>(n\(D- L uD- 2 )) < C* V-ff L~(n) ; ]— ; r 

\n + r 2 J |e log e| 

for sufficiently small e > 0. 

In higher dimensions (n > A), there is a constant C** , independent of e, r*i, r 2 anc 
(ai,a 2; - • • ,a„), satisfying 

(r\T 2 \ 1 
i ~ 
ri + r 2 J e 



W^u\\L^(n\(D 1 uD 2 )) < C**||V-H"||i,«.(n) 
for sufficiently small e > 0. 



r\r 2 \ 1 



ri +r 2 J e 



3 Representation of the potential difference 

Wc introduce a harmonic function h as follows: 



' Ah = 0, inl n \(DiUD 2 ), 

h = 0(\x\ 1 ~ n ), as|x|^oo, 
h\dDi = ki (constant), 

J dD ,d v hdS= fori = 1,2. 

It is essential in this work to construct the function h because of the following lemma. 
Lemma 3.1 For a solution u to ([7]), we have that 



(6) 



M laDi u \dD 2 



[ [d v h)H dS+ f {d v h)H dS. (7) 

JdDi JdD 2 



Proof. With the boundary condition of u on dD\ and dD 2 , Green's identity for H inside 
D\ and D 2 yields 

f d u (u-H)dS=[ d v (u-H)dS = 0, 

JdD x JdD 2 

and thus 

/:= / hd v (u - H) dS + f hd v {u-H) dS = 0. 

JdDt JdD 2 
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Applying again Green's identity outside D\ U D 2 , we have 

= / = (u\ dDi - u| ) - I {d v h)H dS- J {d u h)H dS. 

J dD± J dD2 

□ 

We remark that the above representation ([7]) is observed by Yun [11] for the purpose of 
estimating the stresses between two arbitrary shaped inclusions in R 2 . By constructing a 
harmonic function h and calculating the right hand side of ([7]) , Yun estimated the potential 
difference between two adjacent conductors. 

The idea to establish h is from the basic theory in electrodynamics, and we use several 
times the following property of Apollonius circles. 

3.1 Apollonius Circle in R rt 

For a ball B r (c) in R™ and a point p, |p — c| > r, we have 

r 1 _ 1 

|p-c| |x- R(p)\ |x-p|' 

where R is the reflection with respect to B r (c), i.e., 



for all x £ dB r (c), (8) 



R( P ) = r -^ + c. 
|p-c| 2 

A simple application of Apollonius circle is estimating the potential difference of the 
solution to |T]) for two circles with different radii. 

3.2 Estimates in R 2 

We let 

Di = B ri (d) and D 2 = B r2 (c 2 ), 

where ci = (ri + e, 0) and c 2 = {—r 2 — e, 0), and Ri be the reflection with respect to Di, in 
other words, 

J2i(x) = r f (X ~^ +Ci, i = 1, 2. 

|X Ci I 

Let p-L € Di be the fixed point of Ri o R 2 , then R 2 (p 1 )(~: p 2 ) is the fixed point of 
R 2 o R 1 and i?i(p 2 ) = p x . From ©, 



for x G dD\ , 



l x -Pi| 
|x-p 2 | 



IP2 - ci| 
|Pi - c 2 | 



for xedD 



'2 



2- 



Hence, the solution to (O is 



h ■■= ^- (log |x - Pi I - log |x - p 2 |) = i log ' 



2?r v 1 11 1 Zl/ 2?r V|x-p 2 | 

and, from 0, we have the following proposition. 
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Proposition 3.2 Let H(x\,X2) be an entire harmonic function. The solution u to (QP 
satisfies 

u\dD 1 -u\ dD2 = H(p 1 ) - H(-p 2 ) 

(9) 



= 4S Bl ff(0,0)^£g rv /i + O(e). 
Referring to the mean value theorem, there exists a point x 2 between dD\ and dD 2 such that 



\vu(x 2 )\ > \ dxi H(o,o)\\P^-^=- (io) 

\ n+r 2 Ve 

for any sufficiently small e > 0. Moreover, there is a constat C independent of e, r\ and r 2 
such that 



I|Vu|| l ~(jA(d 1 ud 2 )) < C\\VH\\ L oo ( p l )J~- 



r\r 2 1 



n + r 2 

where ft = -B 4 ( ri+r2 )(0, 0). 
Proof. The fixed points satisfies 



Pl =(2./^^V^+O( e ),0) andp 2 =(-2./-^^+O(e),0 

Therefore, we obtain (|9|) and (fT0|) . By virtue of the argument presented by Bao ct al. in 
[4], the upper bound of the gradient is derived from ([9]). In this paper, the same process as 
this proposition to derive the upper bound of the gradient is also presented in the proof of 
Theorem 12.11 Therefore, please refer to the derivation of the upper bound of the gradient 
in the proof of Theorem 12.41 □ 
We remark that the same gradient estimate has been obtained by Ammari et al. in 
[2] . They represented u by single layer potentials of the Laplacian with potential functions 
defined using Kelvin transform Ri, i = 1,2, and obtained the blow-up rate by investigating 
the potential functions. The novelty of their work is that their estimates is not only for a 
extreme conductivity but also for a finite positive constant. However, as for the extreme 
case, our result provides a much simpler method for obtaining the blow-up rate. 

4 Derivation for Theorem 12.11 

Differently from the two dimensional space where the point charge potential, the logarithm, 
separate the multiplication with ratio p into a sum, we cannot constructed h just with two 
point charge potential functions in higher dimensional space. Therefore, we introduce a 
sequential process to build h. 

4.1 Construct h in R n , n > 3 

Let 

Di = B ri (ci) and D 2 = B r2 (c 2 ), 

where 

ci = (ri + e, 0, . . . , 0) and c 2 = (-r 2 — e, 0, . . . , 0). 



r 



We start from a harmonic function hi } o, defined outside of D\, which is 

1 



/ii,o(x) = 



x-ci|"- 2 ' 



Note that hi t o is constant on dD\. However, it is not constant on dD 2 , and we neutralize 
it by adding auxiliary point charge potential h\.\ to make (h\fi + fei,i) constantly zero on 
dD 2 . From jU, h\,i is defined as 



^i,i( x ) 



'2 



c 2 -ci|/ |x-i? 2 (ci 



\n-2 ' 



where i = 1,2, be the reflection with respect to Z^. On the next step, we add hi f2 to 
(^1,0 + and make (/ii,o + ^1,1 + ^1,2) be constant on d-Di, i.e., 



\ Ji— a / \ 71— z -, 



|c2-ciiy vici-i? 2 (ci)i; ix-Ei^Cd))!"- 2 

Consequently, we construct ft.i )OT , m G N as 



where 



f(fliik)*(ci), ifm = 2fc, fc>0, 

Cl ' m \R 2 (RiR 2 ) k (a), if m = 2fc+ 1, k > 0, 

m 

= n ^> f ° r m £ n ' ( i3 ) 



and 



if j = 0, 



^id-c^r if i = 2 ^^ 1 (14) 

if j = 2k + 1 k > 0. 



|C 2 - Ci, 2 fc| 

Similarly, we define h 2%m , m € N as 



h 2 , m =( g2 , m )- 2 (L-^ , (15) 



where 

f(.R 2 .Ri) fc (c 2 ), ifm = 2ft, fc>0, 

' Ri(R 2 Ri) k (c 2 ), if m = 2fc + 1, fc > 



92,m = J!' 02 ';'' foIm£N > 
J'=0 



and 



>2 



P2,j 



|C 2 - C2,2fc-1 



|Cl - C 2 ,2fe| 



if i = o, 

if j = 2fc, jfe > 1 

if j = 2fc + 1 k > 0. 



(18) 



Since (i?ii? 2 f(ci) G £>i and R 2 {RiR 2 r(c 1 ) G D 2 , we have, for j > 1, 

1 



pi j < max ■ 



.-1,2 r, + 2e 1 + ^ 
By the same way, p 2j < 1 1 1 - ie . Hence, 



where r max = max(ri,r 2 ) 



(&,m 

rn— 



,n-2 



< OO, S = 1, 2 



and the two series ^ =0 )i, jm , i = 1,2, are well defined. To get h satisfying the decaying 
condition, we sum the series /ii, m and h 2>m with different weights as the following lemma. 



Lemma 4.1 The solution to (OJ) is given by 



Mx) = 



i i 



v n-2 



m— 



l x - c i,> 



(2 - n)u> n M 
where iv n is the area of the unit sphere, and 

OO 

Q s =£[(-i) ro (< 



1 71-2 



-Qi£ 



(-l) m (92,™) 



n-27 



771=0 



l x - c 2,7 



1 71-2 



(19) 



,71-2 



771=0 



a = 1,2, 



fe=0 fc=0 

Here, q s , m ' s an d c s m 's are defined by US]), US\), U6]). Ji7| j. 



Proof. Let 



71-2 



771 = 



x - Cl- 



x - c 2 . 



then /i* satisfies (|6]) except the last condition, boundary integral conditions. 
Note that 



1 



d 



1 



(2 - n)uj n J dD . dv |x - p| 



n-2 



e?cr(x) 



0, for p G K" \ Di, 

1, for p G Di, 



(20) 
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and we have 



d v h*(x)do-(x) 



(2 - n)u>„ J 9Dl 

oo oo 



k=0 



k=0 



oo oo oo oo 



fc=0 fc=0 
oo 



QiY2{q2,2k) n + Q2^2{qi,2k+i) 

fc=0 fe=l 

9 I/ ft*(x)dcr(x) 



k=0 k=0 
n-2 



(2 - n)u> n J 9D2 



□ 



Lemma 4.2 Assume that the dimension n is 3 and i/ie distance e is sufficiently small. 
Then, there is a positive constant C independent of r\, r%, and e satisfying the following 
properties: 



Estimates for £™ =0 1s,m 
1 d 



Cd + 
1 1 



C d- 



T |loge| < £>, m < CT— |logc|, 

m— 

T |loge| < f>, m < C^llogel 



(21) 
(22) 



where d = — . 

ri 

Estimates for Q s 



hdTi C dT~r 

1 d „ „ d 
< Q 2 < C 

where Q\ and Qi are defined in Lemma \4- 1\ 

Estimates for Y,T=o( c ^^)(q s ,2k) - J2T=o( c s^+i) (q s ,2k+i) • 



£ < (-l) s+1 



^2(c s ,2k)(q S ,2k) - y^(c s ,2fc+l)(fls,2fc+l) 



Lfc=0 



fe=0 



<CV, 



Now, we are ready to prove Theorem 12. II 
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4.2 Proof of Theorem [2TT1 
Potential Difference 

We first consider the case of H(x\, x%, x$) = x%. Then, Lemma |4. 1 1 implies 



laDi u, \dD 2 



k=0 k=0 

2 ^( c l,2fc+l)(<7l,2fc+l) - -JJ y^(c2,2fc)(g2,2fc)i 



M 



k=0 



k=0 



(23) 



where c s j is the xi-coordinate of c s .j for s = 1,2, j £ N, i.e., 

(c S)j , 0,0) = c S)j . 

Lemma 14.21 allows one to estimate four positive valued terms in the right hand side of (|23[) 
so that 



M = Qi ^ ?l:2fe + Ql ^2 <?2,2/e+l 
k=0 

d 2 



k=0 



1 



(d + i) 2 (d + iy- 

loge|. 



loge 



In total, we obtain 



\dD 1 u '\dD 2 



nr 2 1 
ri+r 2 \ loge | 



In the case of H = a 2 x 2 + 0,3X3, the integration (|7|) is zero; all point charges of h lie on xi 
axis. Therefore, there is no potential difference between inclusions, and we have established 
the estimate for the potential difference between D\ and Z?2- 

Therefore, for H = y —1 ajXj, 



u an ~ u 



sr> 2 



rir 2 1_ 

ri + r 2 I loge| 



Lower bound 

The lower bound is obtained by simply applying the mean value theorem. Since \u\jj 1 — u\d 2 I 
behaves as 1/| loge|, the gradient behaves as l/(e| log e| ) , and, more precisely, there is a point 
Xq between D\ and D 2 satisfying that 



1 



01 



nr 2 



C \n + r 2 J |eloge 



1 



< |Vu(xo)| 



where H = J2t=i a i x i- 
Upper bound 

The upper bound of the gradient is derived by applying the methods presented by Bao et 
al Bi- 



ll 



We assume that r\ > f2, and let 



1 



— A = {x G K 3 | T-i x G Al , i = l, 2. 

n 1 

Note that —A is a unit sphere. Define a bounded domain fi, containing —Di and — A 
independently of e, as the sphere #4(0, 0,0). Now, for a solution u to 0$ for if(x) = 
y?., a^j, define the scaled function u as 

m(x) := — u(nx). 

Then m is also the solution to ([1]) for i?(x) = Y^i=i a i x i with ^ A and ^ A instead of A 
and A- The estimate for the difference of u between 9 A and SA in this theorem yields 



\d(^D 1 ) U \d(^D 2 ) 



1 



Ol 



l + d\\0gS\' 



where 



d = — and 5 



c 

'1 ' 



By the maximum principle, we have 



\u - if||i«.(an) < 



2|Lff| 



L~(S~2) 



< C | 0l 



1 



log<5| 



G2 



«3 



and, as a result, 



02 



03 



To estimate |Vm| on <9(^-A U ^A), we define W3 as in 0]: 



Au 3 = inn\(^AU^A) 
U3 = ^ on d{ ^ A U i A) 
W3 = — u on c3Sl 

We draw the attention of readers to Lemma 14.31 14.41 which are modified from [J] to fit our 
problem. For a reader's convenient, we provide the proofs at the end of this section. 

Lemma 4.3 (J^j) There is a constant C independent of d and e such that 
IIV^ + ^IL-cnu^u^)) < C\ ai \ (^J 

Now, in estimating \Vu\ on <9(^-A U ;^A)- it is remained to be derived an upper 
bound of |Vi>3| on d(— A U — A)- To do that, we define the harmonic function p in 



(24) 



0\(^AU^A) as in 0]: 



Ap = m.Q\{^D 1 \J iflj), 
p = on^ulfla), 
p=l on 9fi. 



12 



Note that V3 = ±||«||.L°°(an)|0 = on d{^D\ U -^D^). Moreover, from the fact that 
V3 = —u on the dfl and the maximum principle, for x 6 Q\(-^-Di U ^j-Z^V we have 
— ||M||L°°(an)P < < ||u||i°°(ar2)p. Therefore, by Hopf's Lemma and the maximum principle, 

ll^ w 3||i«.(n\(JL.i3 lU JL£) a )) < ||u||L-(aO)l|Vp|| i oc( y ( _L DlU _L £ , 2)) . 

We apply the following lemma to calculate || Vu3|| £= o( ^ j_^ lU j_£> 2 ))- 
Lemma 4.4 (J^) There is a constant C such that 

!|Vp|| L oo (a j_ £ , lU9 j_ £ , 2) < c~, 

for e small enough. 

Applying Lemma 14.41 we have 

||V« 3 ||^ (9 (JL Dl uiD 2 )) < C\ (\ ai \ (JL-^ — + | ai | + \a 2 \ + |03|) . (25) 



Two bounds |24]) and (|25j) yield 

||Vu||ioo(a(.D lUD3 )) = ||Vw|| i oo( a( _L £ , lU _L £ , 2)) 

sc (l°-l(TTi)j«^r + >l + N) 
SC "( l<, ' l (^)iTl^ + > l + l ° 3l) 

for sufficiently small e > 0. Since |ViZ| is bounded by |ai| + |«2 1 + la.3 1 , we have 

HV(« - H)\\ L ^ d(DlUD2)) < C (\ ai \ ( J^-) -J-j + i(|a 2 | + loal) 

By the harmonicity of u — H in M 3 \ (£)j U 1)2), ||V(u — -ff)llL°°(R 3 \(,DiU-D.>)) nas the same 
upper bound as the above. The fact of |Vi?| < |ffli| + \d2\ + l a 3| is again used so that 

||Vu|Uoo (R 3\ (Z3lUZ > 2)) < c [\ ai \ {j^r^j + ^l a2 l + H)) ■ 

Therefore, we obtain the upper bound of the gradient estimate. □ 
Proof of Lemma 14.31 

From definition, u + v^ is constant on dj^D\ and ^--Di is a unit ball. By Kelvin transform, 
u + U3 can be extended harmonically to ft \ (— -D1.5' U —D2) where 



— D x s , = {x e — L>i|dist(x,S— D x ) > 5' 
n I ri 1 n 

and 
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Similarly, u + v% can be also extended harmonically to 17 \ (-^-Di^ 1 U -^-D2,5") where 



— D 2 s>> = <xe — D 2 dist(x, d— D 2 ) > 5" 
n I n n 



and 



d + 26' 

Furthermore, by the standard estimate for the extension, there is a constant C such that 

max (u~v 3 )~ min (u - v 3 ) < C\a 1 \ ( 7 ) 1 

n\(iD I(f ,uii) v „) sAf^u^iDy,,) V 1 + d / l lo S<*l 

for sufficiently small e > 0. Note that 

5' w 2<5 and 5" w 25. 

By the gradient estimate for harmonic functions, we have 

d 



l|V(u + «3)||i<x, (n \(-l-13 1 u 1 i.£l 2 )) < Ckil 



+ dj \S\og5\' 



□ 



Proof of Lemma 14.41 

Let pi (z = 1,2) be the solution to 



A Pl = infi\(£A) 
Pi = on a(^-A) 

Pi = 1 on d£l 

Then p = pi on the <9f2 U dDt. The maximum principle yields to pi < p. Note that the radii 
of and -^D 2 are 1 and d respectively. 

Consider the harmonic function u which is the solution to 

Av = inB 4 (0,--- ,0)\B ro (c o ), 
u = on 9B ro (co), 
u = l ondB 4 {0, ■■■ ,0), 

where ro < 1 and |co| < 2. Let the harmonic function w be as 

I (ox=s k=s) 1 ( 1 W phr! for 71 > 3 

{ (log 2 - log r ) (log |x - c | - log r ) for n = 2 

Then, 

ic>l = Don as 4 (o, • • • ,0) 

and 

w — v on dB ro (c ). 
Then, there is a constant C, independent of ro and Co, satisfying 

l|Vw|| L oc (ai 3 ro(co)) < C-. 
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By Hopf 's Lemma, we have 

—d v w > —d v v > on dB ro (c ), 
and there is a constant C, independent of Co and ro, satisfying 

-d v w < C-. 

Therefore, we obtain 

\\Vv\\L~( a B T0 ( C0 )) = -d v w<C±-. (26) 



It follows from (|26|) and Hopf 's Lemma that 

IIVpHioo^jLoj) < llVpiHioo^jL^) <C 



and 



□ 



4.3 Proof of Lemma 14.21 

We begin by considering the position sequences Ci jm and c 2jm , because the quantities like 
qi t m and pij used in the derivation are yielded by these position sequences c,. m . Referring to 
the relation (p~2]) . the successor Ci jm+2 to ci. m is determined by the twice refection i?i o R 2 . 
For x = (x, 0, . . . , 0) G Di, the twice reflected point i?i(i? 2 (x)) = (x', 0, . . . , 0) is given by 

r 2 

x = n + e ■ 



n + r 2 + 2e 



X+T2+6 

For the sake of convenience, we assume that 

e j r 2 _„j .. c i,j 



6= — . d= — , and y, ,= j € N, (27) 

n n ri 



where Ci ,• is the Xi-coordinate of Ci Then, we have the relation 



d+ (1 + M)5 + 25 2 ni d+(3 + d)5 + 2S 2 

d + 25 

AdS + i{l + d)8 2 + 25 z 



V2kV2k - 2 + l + d + 28 V2k l + d + 25 V2k ~ 2 



I + d + 25 



= 0. (28) 



Let p = (p, 0, . . . , 0) € Di be the fixed point of (Ri o i? 2 ), then ^- is the limit point of y 2 fe 
and satisfies 

/_p\ 2 2{d-l)5 fp\ 4dS + 3(1 + d)5 2 + 2(5 3 _ 
V^J + l + rf+2(5 VnJ 1 + d + 2(5 ~ ' 

and, as a results, 

^= 2 y^ + 0(*). (29) 
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Here, we have a constant C independent of d so that 

\0{y/5)\ < CV5 for sufficiently small 5 > 0. 



It follows from (121) that 



( p 

V2k 

V n 



p\ ( d + (1 + 3d)5 + 25 2 p\( p 
V2k~2 + : , , , of 1 2/2fc 



n / V ri / V 1 + d + 25 ?*i / V n 

d+ (3 + <)£ + 2<5 2 p\/ p 



For simplicity, let 



l + d + 26 ri )\ V2k 2 n 
P 



^2fc = 2/2fc 



then 

1 ^rf + (1 + 3rf)(5 + 2<5 2 p^ 1 + (3 + d)<5 + 2(5 2 p^ 1 

Further, let 



V 1 + d + 2<5 n/ Z2k-2 ^ 1 + d + 2<5 ri/ z 2 k 



( d+(l+3rf)<5+2<5 2 \ _p_ 
. I l+d+25 J + n 

A = —, r and B 



( d+(3+d)S+2S 2 \ 


_ JL 


\ l+d+25 J 


ri 



o ( (d-l)S , p_\ 

* \ i+d+25 T n y 



this can be rewritten as 

/ 1 \ / 1 



— + 5 = A + B 

\z 2 k ) VZ2fe-2 



Therefore, the sequence y 2 k is expressed as 

U2k = z 2 k H 

n 

i 



(i + B)A fe - b n 

- I + JL 

(t+ct + B ) Ak - B 



(30) 



4.3.1 Estimates for yj 

We simplify ((3"0)) under the assumption that S is sufficiently small. Our strategy is to choose 
an appropriate N so that ^2 k<N is dominant in the following series calculation, for example 
in XX<7i.2m — <?i,2m+i)i and estimate the series separately for smaller and larger sub-indices. 

From now on, we use the big O notation frequently. The equation f = g + 0(d) means 
that there exist a constant C independent of 8 such that |/ — g\ < CS for small enough 
S > 0. In this paper, C is assumed additionally to be independent of d and k as well. We 
define 0(V~5) similarly. 

The expression (|3U)) of yj is too complicated to well describe the dependency of yj with 
respect to d and S. Thus, a simplified expressions is established, provided that the distance 
5 is small enough. As for our strategy, we choose an appropriate number N ~ y/e so that 
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the sequence terms of k < TV are dominant in the sequence yk , and thus estimate j/2fe hi two 
cases of k < TV and fc > TV separately. 

From the definition of A and B, we have 



.4 



d 
1+d 



0{5) + "=l + 2 d -±±^ + d 4±0(6) 



d ri 



, d + l R d+l^.„. 
1 + 4a/ — — VS + ——0(6), 



and 



y/SB 



1 d+l 



d+l 



o(Vs). 



By a standard argument, one can show that for x £ (0, 2) and (1 + x) k < 2, 

(l + x) k < l + kx + k 2 x 2 . 

log 2 



For k < ^a/^^j, we have (l + 8^^Vsj < 2. This yields 



A k = 1 + k 4 



d+l rz d+l_.„ 
— j— V (5 + —j—O (J) 



d+l 



/fc 2 0((5) 



Hence, the estimate implies 

V2k 



(B + 1 + 0{VS)) ( 1 + AkiJ^VS + ^k 2 0(8) J - S 



l + fe i±i + (d±l)^2 0l( ^) n 

Here is a constant C\ > independent of (5, d and A; such that 

\Oi(VS)\ < dVs. 

We take the integer TV as 



TV = min < n > 

nGN 



1 log 2 / d 1 



d 8 Vrf+1^ 



Then, for k < TV, 



fc(cZ + 1) + d 



(31) 
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and 

1 



2/2fc+i = — - r 2 - e 
n V 

= -d-5+ 



ri\ r x y 2k + r 2 + e 

d 2 



k(d+i)+d T » T V 3TT (^) 
lk + 1 M d+1 y ' 



i. - (32) 

In the case of k > N, we use the fact that the sequence j/i,2fc is decreasing to i.e., 

2/0 > 2/2 > • • • > 2/2fc > 2/2(fe+i) >•••>£• Here ((311) yields yi. 2N = \f^T°(^)- From 
the estimate (|29p for we have 

C^S6>y»>^r S (33) 

and 

Therefore, all estimates (|3Tj) . (|32|) . (|33|) and (|34[) for y 2k and j/2fc+i are obtained. 

4.3.2 Estimates for J2m=o 9s,m 
We consider the estimate for 

rn 

qi.;n = Y[ /'' ,'• 

For fc < N, from (JIJ), (glj and (32]), we have 

^ = jtst^ = (fcTixW 1 + Twry ^^' (35) 

pl - 2k+2 = i + s- V2k+1 = + + + + V ^TT 0( ^ }) ' 

and 



In the case of m < TV, they lead to 



m—l 



=( II Pl^ + lPl,2 fe+2 ) = m(d+1) + d (l + V ~ O^)" ^ 



fc=0 

d 



d+l 



m{d+ l) + d 

?l,2ro+l = (<Zl,2ro)(Pl,2ro+l) 

d / ld+1 



(m + l)(d+l) 



) 

(i+ m y?no(V3)) 



18 



Thus, there exists a positive constant C such that, for s = 0, 1, 

1 d 1 „ d 1 

< 5l,2m+s < C- — , 1 < m < N 



Cd+lm+1 ~ ~ d+lm+1 

and, therefore, 



-| log <S| < Yl ( 

9l,2m + <7l,2m+l) < d~+T ^ ' 



In the case of m > iV, (O and (O yield, for s = 1, 2, 



C 

Then it leads to 



<p s ,m<C 1 , m>iV. 



1+ / minaMrf) ^ " 



/ , (91, 2m + <Zl,2m+l) — \Ql,2N + ?l,2JV+l) / ,l 
— (qi.2N + Ql,2N+l) 



d 1 1 



d + 1 N / min(l/dj) ^ 

~max(d,l). (37) 
Therefore, we obtain (f2Tj) . and replacing c? by 3 we also have 



4.3.3 Estimates for Q s 

We consider Qi = X)m=o (#i,2m — <7i,2m+i) • From definition, gi, m has the decreasing property 
as 

9l,2m+l = (pi,2m+l)(<Zl,2m) < 91,2m) 
<Zl,2m = (/>l,2m)(5l,2m-l) < <?l,2ro-li 

and therefore 

< H (91,2m - 3l,2m+l) < gi,2JV < C ( d+1 )(jV + 1 ) - ( 38 ) 

This means that X^Tn>A r ('?i,2m — <7i,2m+i) shrinks to as 6 goes to 0. 
On the other hand, it follows from (|3l)j) and (|3l)|) that 



0<m<N m<N 



(9l,2m — <Zl,2m+l) — J]] 9l,2m (1 — Pl,2m+l) 

m<JV 



m<N 
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Taking an advantage of a strict decreasing sequence, we get 
1 / 1 



&tt+«)< 51 



< 



where 



0<m<AT 



M 



(m + l)(d + l)(m(d + 1) + d) d+1 



21 



o (t + l)(d + l)(i(d+l) + d) 



d+1 
d 



d+1 



t+1 



dt 



d+1 



log^-+log(^ rTr . 



As has been mentioned, S is assumed to be small enough so that 

^- log (^vrr^2 log ^r- 



It leads to 



1 d d+1 I, ^ 3d,, d+1 1, 

'log— ^ + 3) <gi<^^^(io g - 



4d+ 1 



Note that 

Therefore, we obtain 



d d 1 ~ ^ 2d+l 
log(l + ar) < x, for x > 0. 



< Qi < Ci 



Ci d+1 



d+1 



Similarly, replacing d by 4, we also have 



Id „ d 

< Q2 < C 2 



C 2 d+1 



d+1 



4.3.4 Estimates for Efelo( c ^2fe) (<? s ,2fe) - Lfclo( c ^2fc+i) (g s ,2fc+i) 
The lower and upper bounds of 



^]( c s,2m)(9s,2m) — (c s ,2ro+l) (gs,2m+l); S — 1, 2 
m— m— 

are established here. From (|2"7|) . (f5T|> . I|32p and (|37p. we calculate 

— ^ (Cl,2m) (gl,2m) = X] (#l,2m) (gl,2m) + E (i/l,2m) (gl, 

d 



2m 



m=0 



0<m<Ar 



0<OT<JV 

< C. 



m>N 
2 /d+1 



m(d+ 1) + d 



21, 



0(V5) £ gi, 



m>N 
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Moreover, 



Similarly, we have 



— YYci, 2 fe)(9i,2fe) > (yi,o)(9i,o) = 1- 



k=0 



Y oo 



fe=0 



and, in total, we have 



1 1 

C ~ n 



^(ci,2fc)(?l,2fc) ~ 53(ci,2fc+l)(9l,2fc+l] 



.fc=0 



fc=0 



By the same way, 



1 1 

c ~ ~V 2 



^2{c2 : 2k){q2,2k) - y^(C2,2fc+l)(g2,2fc+l; 



< c. 



< c. 



k=0 



k=0 



□ 



4.4 The Derivation for Theorem 12.31 

Lemma 4.5 Assume that the dimension n > 4 and i/ie distance e is sufficiently small. 
Then, there is a positive constant C independent of r\, r%, and e satisfying the following 
properties: 

• Estimates for £m=o<m ; 



m— 



m=0 



where d = — . 
Estimates for Q s 



Qi ^ , 1 , and Q 2 " 



d+1 



d+1 



where Q± and Q2 are defined in Lemma \4- 1\ 

Estimates for £^L (c s , 2 fc) (g s ,2fc) - £^L ( c ^2fc+i) (?s,2fc+i} 



(-1)- 



-1 



^(C s ,2fe)(<7s,2fc)" 2 ~ y^(c 8 ,2fc+l)(gs,2fc+l)' 



,fc=0 



fe=0 



r s for s = 1, 2. 



Proof. Let JV be as chosen in the proof of Lemma 14.21 We have shown in Estimates for 
£m=o Qs,m that for j = 1, 2, 



9l,2m+j 



d 1 
d + 1 m + 1 



for 1< m < N 
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and qifi = 1. This yields 

V <7"- 2 ~l 
0<m<JV 

and by the argument of (|5T)) . we have 



i(n-a) 



— ^l,2iV + «1,2AT+1 



d l x "" 2 



d+liV 



Since N increase at the rate of as S goes to zero, Y^ m >N{li,2m + 9i.2m+i) shrinks to 
zero as 5 goes to zero. Therefore, we obtain X)m=o ?"m — 1; an d similarly X)m=o 92 m 2 — 1- 
Now we consider Q s = J2m=o(^^-) m ( ( ls,m) n ~ 2 for s = 1, 2. We note that the tail sum 



^ (-lrfa,*,)"- 3 ! < i e <2»+<2™ + ii 



m>2N m>N 



and the upper bound in the right hand side has been shown above to shrink to zero as 8 
goes to zero. In this respect, we estimate only S to <2jv-i( — l) m (9i,m) n-2 f° r Qi- From 
definition, we have 



£ (-l)-(g 1>ro )«- 2 = J2 (<Zl,2m)"- 2 (1 - (P2, 2m+ l)- 2 ) 

( ^TTTTTZT/ ) (i-^m+i) 



m<2iV-l m<JV 

n-2 



m<N 



m<N 



d 


m(d - 


f l) + d 




d 


?7i(c? - 


f l) + d 




d 



^(wrf+-n + J ( ( m + l)(d+l) + ^(d+l) ^ 



rn<AT 

n— 2 



m(d+l) + dy \(m+l)(d+l) 



1 d n ~ 2 
d+l + 2-" (d+ l) n - 1 (m + l)™" 2 

Km<JV \ / \ / 

1 



d + l 

Therefore, we obtain Qi ~ ^^i> and replacing d by ^, also have Q2 — jfpj- 
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We consider the last estimate. By Lemma T4. 21 we have 

i < ^V- 1 



< 



< 



1 

i 

n 



,k=0 k=Q 



^2(ci,2k)(qi,2k) - y^(ci,2fc+l)(gl,2fc+l) 



,k=0 



k=0 



Therefore, we have 



^(ci,2fc)(gi,2fe)" 2 - ^2(ci,2k+i)(qi,2k+i)' 



Lfc=0 



fc=0 



and similarly 



^2{c 2 ,2k)(q2,2k) n - y^(c2,2fc+l)(g2,2fc+l 



, n-2 



.fc=0 



fc=0 



'2 



□ 



4.4.1 The proof of Theorem! 

We first consider the case of H{x\,X2 1 - • ■ ,x n ) = x±. Then, Lemma |4 . 1 1 implies 

Q 



u \dDi u \dD 2 m 



= -T7 y^(ci,2fc)(gl,2fc)" + "77 y^(c 2 ,2fc+l)(g2,2fc+l 



M 

fc=0 fc=0 



Q2 2 Ql 

-JJ^2(ci.2k+l){qi,2k+l) n - JJ^2(c2,2k)(q2. 



2k) n ~ 2 - (39) 



fc=0 fc=0 

where c Sj j is the xi-coordinate of c S J for s = 1, 2, j e N, i.e., 

(c s ,j, 0, • ■ • ,0) = c„ t j. 

Lemma 14.51 allows one to estimate four positive valued terms in the right hand side of 
so that 



M = Q2 2J 3l,2fc + Ql X] ?2,2fe+l 



A;=0 



fe=0 



In total, we obtain 



U \dD! U \dD, 



ri + r 2 
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In the case of H = ^™_ 2 a i x i> the integration (J7J) is zero, because all point charges of h lie 
on x\ axis. Thus, there is no potential difference between inclusions. Therefore, we have 
established the estimate for the potential difference between D\ and D 2 . 

Now, we consider the upper bound of |Vu| when -ff (x) = Yh=i a i x i- To do so, we pursue 
the argument similar to Theorem 1 2. II We thus assume that r% >r 2 . Let ^r-Di and j^D 2 be 
the same as defined in the proof of Theorem 12. 11 We choose the sphere B^O, . . . , 0) as the 
domain Q containing D\ and D 2 independently of any small distance e > 0. We consider 

m(x) := — u(rix). 
ri 

Then, u is the solution to ([1]) for -ff (x) = Y2i=i a i x % when the inclusions are ±D X and ±D 2 
instead of D\ and D 2 . It follows from the estimate for the difference of u between D\ and 
D 2 that 



«i 



1 + 



where d = £a and S = ~. As defined to (0.9) in [I], «3 is defined as follows: 



Av 3 = inn\(i J D 1 U^ 2 ) 
w 3 = on U ±D 2 ) 

1)3 = —u on dfl 



Lemma 14.31 and 14.41 have been used in the proof of Theorem 12.11 By the help of them, we 
can obtain 

l|V(w + w 3 )|| i ^ A( j_ DiU _L D2)) < C\a x \ f — 



and 

\\Vv 3 \\ L ^ {d{ j_ DiU j_ D2)) < C- d Mai| (j 

Two bounds lead to 

||Vu||£«(a(x) 1 ui3 2 )) = |jVu|| L oo (a( j_ £)lU j_ £l2)) 





for any sufficiently small e > 0. Since |Vff| is bounded by ^ZILi l°t|) the harmonicity of 
u - # in W l \ Q. leads to 

||Vu|Uoo (R n\ (i3lUD2)) < C ^|0l| ^-^-^ i + I^loiM . 

Therefore, we obtain the upper bound of the gradient estimate. In addition, the lower bound 
of the gradient is immediately derived from the mean value theorem. □ 
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4.5 The proof for Theorem 12.41 

Without the loss of generality, we may assume that H (0, • • • , 0) = 0. Then, Lemma 14.11 
implies 



\6D X U \dD 2 



yrc-2 Ql T w w \n-2 



k=0 



M 



<l|VF|| L =o ( n)|^^(ci, 2fc )( g i 



,2fc 



v n-2 Q 

M 



7 y^( c 2,2A;+l)(g2,2fc+l 



,n-2 



fc=0 fc=0 
oo 



^ y^( c l,2fc+l)(gl,2fc+l) - J7 y^(c2.2fc)(g2,2fc) 



fe=0 



fe=0 



Owing to Lemma 14.21 and 14.51 we can choose a constant C independent of e, r\, r 2 and 
|jV J ff|| L = ( n) satisfying 



I an. 



<C||VfT||io 



^1^2 



tt^ — r if n = 3 

I loge| 

1 if n > 4. 



Here, we would establish the upper bound of the gradient. The method is the same 
as presented in the proofs of Theorem 12.11 and 12.31 We thus assume that r\ > r 2 . Let 
— Di (i = 1, 2), u, d and S be as defined in Theorem 12.11 and 1 2.31 Then u is the solution to 
CL]) for 

H :=- H(nx) 
r i 

instead of H, replacing Di by ^--Di for i = 1,2. We consider := {x 6 1" | rix G f2}, 
because fl is the ball B i ^ ri+r ^(0, ' ' ' >0) m this theorem. So, the radius of the ball ^fi is 
between 4 and 8. Note that 



||VF|U~ ( n) = IIVFILoo^). 

We thus choose the harmonic function W3 satisfying the same definition as Theorem 12.11 and 
2.31 replacing fl by ^Q. Similarly to Theorem 12. II and 12.31 owing to Lemma r4.3b nd 14^41 or 
the lemmas presented in [J], we obtain 



\\V(u + v 3 )\\ Lao{ j_ aV j_ DlU j_ D2)) < C||Vff|U~ (n) — • i f 



if n = 3 
if n > 4 



and 



IIv^il-^jl^u-^)) < c-\\vh\\ l ^ u) y—d+i 

Two bounds above yield 



||VM|| L oo (a( _L £ , lU J_ £ , 2)) < C\\VH\\ L o, {n) I - 



m — it if n = 3 

l<5 logA| 

i if n > 4 
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By the harmonicity of u - H in W l \ (^-fi), we have 

|| V(U - H")|U=o(R«\(D 1 uB 2 )) = II _ - ff )!lL~(R"\(J r _D 1 uJ r n 2 )) 

Furthermore, this leads to 



|Vu|| L oc (n \(.D 1 u.D 2 )) < C\\^JH\\ L ^ m 



r i r 2 \ f ifn = 3 



Ji+r 2 ) \ \ ifn>4 
Therefore, we have completed the proof. □ 
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